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Capillary forces exerted by liquid drops caught between crossed cylinders: 
A 3-D meniscus problem with free contact line 
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Abstract 

The Young-Laplace equation is solved for three-dimensional menisci between crossed 
cylinders, with either the contact line fixed or the contact angle prescribed, by means of 
the Galerkin/f inite element method. Shapes are computed, and with them the practically 
important quantities: drop volume, wetted area, capillary pressure force, surface tension 

force, and the total force exerted by the drop on each cylinder. 

1. Problem statement 


A liquid drop between solid fibers {Fig. 1) constitutes a three-phase system whose 
thermodynamic aesuription may become quite involved (cf . Huh 1969) . However, even without 
gravity, if the fluid-solid contact line differs from a circle, or two coaxial circles, 
the drop surface generally forma a three-dimensional meniscus of constant curvature, dif- 
ficult to approximate from experiment and non-trivial to compute from theory. In large 
drops, gravity joins surface tension in the molding of shape and the situation becomes 
even more complex. The surface energy of the system in the absence of gravity is 
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Here the solid surface is rigid and the solid-liquid vapor interactions are described by 
the constant interfacial tensions c T , r , a T _, and 

LV Eb bV 

Equilibrium drop shapes make the surface energy stationary, i.e. 6E * 0 for all pertur- 
bations that leave the drop volume fixed, 6V * 0 , and either the contact line fixed, 

6x L sv * Q, or the contact angle unchanged, ( $ L V ‘ Os * cos 6 t where 6 is the contact 
angle) . These conditions define a variational problem equivalent to the Young-Laplace 
equation of capillarity — a normal stress balance — with appropriate boundary conditions 
at the contact line: 
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Eq. [2] makes plain that the vanishing variation forces the Young-Laplace residual, 

2Hj y -(py - p L )/0 LV , and likewise the Young-Dupr4 residual, cosB-(osv * °LS)/°a,V • to 
orthogonal (as viewed in the appropriate function space) to admissible ncrmal perturba- 
tions . 
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A drop of liquid trapped between two perpendicular cylinders of equal radii is of in- 
terest. When the cylinders touch, their point of contact is taken as the origin (Fig. 2). 
Spherical coordinates are used to obtain a single-valued representation of the liquid- 
vapor interface. A sphere of radius R, equal to the cylinder diameter, is chosen as the 
case surface and is centered at the origin. 

The unknown position of the interface, measured in the units of R, is 

* LV - f (e ,♦) e r [3] 


Analogously, the interface perturbation is 
6x lv ■ C(8,*> e r 


[ 4 ] 


Furthermore, in the case of spherical representation, Weincarten’s formulas relate the 
mean curvature of the interface to the divergence of the unit normal on the unit sphere ft 
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This, together with the divergence theorem, leads to the weak form of Eq. [2] : 
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2. Galerkin/f inite element method 

The main steps of the finite element algorithm can be outlined as follows: 

(i) The problem is nonlinear and so an initial estimate of drop shape ( and contact 
line poeition) is needed that falls within the domain of convergence of the 
iteration method. Estimate from either analysis of a limiting case or experi- 
mental observations or a blend of both. The last leads to the inside surface of 
a email toruc as an initial guoss. 

(ii) To discretize the problem, partition the corresponding spnerical domain (Fig. 3) 
into curvilinear quadrilaterals between the equally spaced spines $ * constant 
(cf. Kistler 1981). The sp .nes remain fixed but the nodes of the quadrilaterals 
car. move along them. Let the number of nodes in the partition be N. 

(iii) Construct a finite element basis function f 1 (6,$) for the subdomain around each 
node (cf . Brown, Orr, and Scriven 1979) : 

- choose the biquadratic polynomial on the (£,n)-unit square, 

- map each quadrilateral isoparametrically (8,$) (£,n) onto the unit square. 

This procedure transforms the original, free boundary domain ft in (6,$)- 
coordinates into a fixed square in the map (£,n) (Strang and Fix 1973). 

(iv) Approximate the drop shape as 

f '9,*) -Ea i f i j%(£,n),$K,n)J 

- nodal values a i are the coefficients to be found. 
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(v) Approximate the shape perturbation as 

N if 

5(0,$) = z Y i ¥ [^0(5, n) , <M5,n) 

take as arbitrary coefficients. Because the y- 's are arbitrary, the 

Galerkin weighted residuals — to which Eq. [6] transforms under (iv) and 
(v) — must vanish at equilibrium. This brings out the direct link between 
the Galerkin and var national approach. 

(vi) Solve the resulting N nonlinear algebraic equations for the coefficients 

using Newton's method. 1 

update the original domain in each iteration 
- use the Jacobian for continuation in the parameter (py - Pi,)/ a LV or v * 

(vii) Terminate iteration when the L^ norm of residuals, i.e. the largest residual, 
is smaller than a preset value, e.g. 10”^ (as was actually used). 

The algorithm was programmed in Fortran and executed on a CDC CYBER-74 computer. U 
took ca. 2 sec/iteration for 169 unknowns and three to at most five iterations to converge. 

A sequence of drop shapes in order of increasing volume is shown in Fig. 4 
(R(p v - P L )/0Ty = LAMBDA) . As capillary pressure decreases, the liquid forms drops of 
increasing volume and finally encircles both cylinders. The plotted drops are stable, 
including the largest one — the surface of which covers more than half of the cylinder 
cross-section. 


3. Practical quantities 


Drops of liquid caught between crossed fibers of non-woven fabrics draw them together 
by capillary action and, when the liquid solidifies, fasten the fibers together. Similar 
phenomena occur in paper, where water at times forms droplets by capillary condensation 
from humid air. Liquid drops between fibers are also present at certain stages of oily 
soil removal by detergent action. Furthermore, drop behavior between crossed cylinders 
is important to direct measurement of adhesive forces between solid surfaces in the pre- 
sence of capillary-condensed liquids (Fishe- and Israelachvili 1981) . These and other 
applications call for computation of the following practically important quantities: 
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their sum: total force holding crossed cylinders together. 

Fig. 5 shows that the total force on the cylinders increases with the wettability of 
the solid by the liquid. As can be seen from Fig. 6, the components of total force vary 
widely with the drop volume but their sum increases slowly and approaches an asymptotic 
limit as the volume shrinks to zero. That limit is 4 tt a^y (R/2) cos 6 (Fisher and 
Israelachvili 1981), and becomes a good approximation for R(py * PL^/°LV % 10 4 , i.e. for 
very small drops. If the mean curvature is taken as the independent parameter, the drop 
volume grows extremely sensitive to it as the contact angle decreases (Fig. 7), whereas 
the total force again approaches the limit (Fig. 8). 


Summary and conclusions 

What is described here is an extension cf earlier analyses of three-dimensional menis- 
cus shapes (Orr, Brown and Scriven 1977) . 

Besides yielding equilibrium shapes of drops with fixed contact linee s as in past 3-D 
analyses, the algorithm used here also gives shapes of drops with a prescribed contact 
angle , which is sometimes a closer approximation to reality. 

Computations for other cases of 3-D menisci making a prescribed contact angle can be 
treated similarly. Augmented by a block-Lanczos method for solving the related 
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eigenproblem (see Brown and Scriven 19BCH , the algorithm tests stability with respect to 
admissible perturbations with the contact line fixed. The case of prescribed contact 
angle is more difficult but is also being treated in the continuation of this work. 

The results show that total capillary force between cylinders increases with decreasing 
contact angle, i.t* with better wetting. Capillary force also increases with decreasing 
drop volume, approaching an asymptotic limit. However, the wetted area on each cylinder 
decreases with decreasing drop volume, which raises the question of the optimum drop voi~* 
ume to strive for, when permanent bonding is sought from solidified liquid. For then the 
strength of the bond is likely to depend upon the area of contact, which is the wetted 
area when the bonding agent was introduced in liquid form. 
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Fig, 1. Liquid drop caught between crossed cylinders 
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Two cylinders of equal radii crossed 
at right angle — coordinate system 
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Fig . 3 

Typical computational domain in spherical coordinates 
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Fig. 4 

Sequence of drop shapes in order of increasing volume. 
Dimensionless capillary pressure R(Py - P^)/o^y * LAMBDA. 
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Total force exerted by drop Components of the total force 

on each cylinder vs. drop voJume. vs. drop volume. 
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Fig. 8 


Drop volume vs. 
capillary pressure 


Forces exerted by 
drop on each cylinder 
vs. capillary pressure 
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Fig. 9 

Wetted area on each cylinder vs. drop volume 


